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Abstract. In this contribution, we present the StaggerGrid, a collaborative project for the 
construction of a comprehensive grid of time-dependent, three-dimensional (3D), hydrodynamic 
model atmospheres of solar- and late-type stars with different effective temperatures, surface 
gravities, and chemical compositions. We illustrate the main characteristics of these 3D models 
and their effects on the predicted strengths, wavelength-shifts, and shapes of spectral lines, 
highlighting the differences with respect to calculations based on classical, one-dimensional, 
hydrostatic models, and discuss some of their possible applications to elemental abundance 
analysis of stellar spectra in the context of large observational surveys. 



1. Introduction 

The Gaia mission [HE] will measure high-precision parallaxes and proper motions for 10^ galactic 
stars down to apparent magnitude V = 20, as well as radial velocities for 10* stars with V<13. 
In connection with the mission, a number of ground-based surveys are being planned that will 
carry out observations to complement the Gaia data. Over the next five years, the upcoming 
Gaia-ESO public survey, for instance, will acquire high-resolution spectra for more than 10^ 
galactic stars using the VLT/FLAMES multi- fibre spectrograph. The goal of the survey is to 
homogeneously derive chemical abundances for all these stars; combined with astrometry from 
Gaia, this information will allow to trace the most detailed and extensive chemo-dynamical map 
of the stellar components of the Milky Way. 

Processing the enormous amount of data from this and similar surveys and extracting from 
them stellar elemental abundances will require fast access to grids of model stellar atmospheres 
and synthetic spectra covering the relevant range in terms of stellar parameters and chemical 
compositions. Grids of traditional one-dimensional (ID), stationary, hydrostatic model stellar 
atmospheres such as MARCS [3l H] or ATLAS O [6] are already available or are being re- 
computed together with atlases of synthetic spectra to match the specific needs of the various 
surveys. The main advantage of classical model atmospheres is that the simplifying assumption 
of a ID stratification allows to invest the computational resources on the solution of the radiative 
transfer equation for a very large number (typically 10^ or more) of wavelength points. However, 



ID models can only treat convective energy transport in an approximate manner via recipes 
such as the mixing-length theory (MLT) [7| or the full-spectrum-of-turbulence (FTS) model 
[S], which are all dependent on a number of free parameters. In late-type stars, the convective 
flows reach the stellar surface, affecting the atmospheric layers and, consequently, the actual 
spectral energy distribution in the emergent radiative flux. Furthermore, bulk gas flows in 
stellar atmospheres and associated Doppler shifts also affect the broadening, shapes, wavelength 
shifts, and strengths of spectral lines. It is therefore important to properly account for their 
effects in order to extract accurate and precise elemental abundances from the analysis of stellar 
spectra. Because traditional stationary, ID, hydrostatic model atmospheres of late- type stars 
lack a consistent description of atmospheric velocity fields, Doppler broadening of spectral lines 
is modelled in ID analyses by introducing additional free parameters such as micro-turbulence 
and macro-turbulence that generally need adjusting and tuning on an individual star basis. 

In more recent years, on the other hand, a lot of efforts have been invested in the development 
of realistic, time-dependent, three-dimensional (3D), radiation-hydrodynamic simulations of 
stellar surface convection that can be directly applied as 3D model stellar atmospheres for 
spectral synthesis calculations (see [9] for a recent review). Three-dimensional simulations of 
stellar surface convection successfully reproduce important observational constraints such as the 
spatial properties and temporal evolution of the solar granulation pattern [IHl Ell H] j centre-to- 
limb intensity variations at the surface of the Sun \13\ [T3] as well as other nearby stars \15\ [16] , 
as well as the detailed shapes and wavelength shifts of spectral lines in solar- and late-type stars 
[13 EHl IS! • At present, a number of codes suitable for constructing realistic 3D hydrodynamic 
model atmospheres are available and actively developed. A few prominent examples in this 
respect are the Stagger [20J, Bifrost [21], CO^BOLD [22], and MURaM [23] codes. Three- 
dimensional models have recently started to be employed for spectroscopic stellar abundance 
analyses [211 [25l [26l |271 |28]. Similarly as in the case of ID models, organized grids of 3D 
model atmospheres are currently being developed and becoming available for this purpose. The 
recently presented CIFIST collection of CO^BOLD model atmospheres [29], for instance, is 
a pioneering example of such a grid. In this contribution, we will present the general outlines 
of the StaggerGrid, an alternative project for the construction of a grid of 3D model stellar 
atmospheres of late-type stars with the Stagger-Code, and we will discuss some of its possible 
applications in the context of spectral line formation and elemental abundance analysis. 

2. The StaggerGrid 

We are using a custom version of the MPI-parallel Stagger-Code to carry out time-dependent, 
three-dimensional, radiation-hydrodynamic simulations of convection at the surface of late-type 
stars for a range of effective temperatures, surface gravities, and metallicities. The general 
outline of the grid in stellar parameter space is shown in Fig. [TJ 

The simulations are of the so-called star-in-the-box kind: each simulation's domain is 
a rectangular, three-dimensional volume located at the surface of the star, with periodic 
boundaries horizontally and open boundaries vertically. The mass, momentum, and energy 
conservation equations for a compressible, viscous flow are discretized and solved on a Cartesian 
numerical mesh with 240x240x240 grid-points as default numerical resolution. The code uses 
artificial viscosity operators to deal with the effects of numerical diffusion. The parameters 
for these operators have been tuned using solar surface convection simulations to minimize 
such effects without causing undesired or excessive smoothing of shocks. Once tuned, these 
parameters are kept fixed for all models. Changing these parameters does affect the sharpness of, 
e.g., temperature and density inhomogeneities but essentially does not alter the main features of 
the simulations, such as the basic morphology of the convective flows or the average temperature 
and density stratifications. 

We choose the physical size of the domains so that the simulations typically host about ten 
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Figure 1. Overview of the StaggerGrid: the dots indicate the effective temperatures (Teg) 
and surface gravities (log g) of the 3D surface convection simulations being computed for the 
project. The coloured circles represent the different compositions (scaled solar metallicity |14) 
with [Fe/H] ranging from +0 to —3, with a a-enhancement value of [a/Fe]= +0.4 for the metal- 
poor models). The empty squares indicate the T^s- and log (7- values adopted for the surface 
convection simulations of the Sun and other standard stars important for stellar spectroscopy 
(see also Tabled]). 



convective granules at any given time and cover about twelve pressure scale heights vertically, 
extending from logTR,oss<— 4 to logrRoss>+6 in terms of Rosseland optical depth. In the setup 
currently used for the StaggerGrid models, we assume a constant vertical gravitational 
acceleration in the simulation box and neglect sphericity effects. The physical size of the 
computational domain of low-surface-gravity StaggerGrid models is, however, still sufficiently 
small compared with the stellar radius, implying that these approximations are acceptable. At 
the bottom boundary, located deep below the surface, in the convectively unstable layers, we 
impose constant pressure across the whole layer and require the inflowing gas to have constant 
entropy per unit mass. 

Our goal is to construct model atmospheres that can be used for accurate abundance 
determinations, so we try to include as realistic input physics as possible in the simulations: 
we implement a state-of-the-art equation-of-state [30] and up-to-date continuous [3] (see also 
Trampedach et al., in prep.) and line opacities 

In order to model the temperature stratification in the outer stellar layers properly, it is 
paramount to account for the energy exchange between gas and radiation. At each time-step 
during the simulation, we solve the radiative transfer equation along the vertical as well as 
eight other inclined directions (two 6- and four (/>-angles) using a Feautrier-like method |31j 
and compute the necessary radiative heating rates for the energy conservation equation. The 
choice of a Feautrier-like method implies that each direction is effectively counted twice, once 
for outgoing rays and once for incoming rays. 



Temperature gradients at the surface of red giant models can become very steep. In general, 
a fixed geometrical depth scale does not allow to properly resolve these gradients, which can 
give rise to numerical artifacts in the outgoing intensity pattern at the surface. In our version of 
the Stagger-Code, we have therefore implemented an adaptive depth scale to solve this issue 
(Collet et al., in prep.). 

In order to reduce the computational burden for the solution of the radiative transfer 
equation, we approximate the source function with the Planck function at the local temperature 
{Su = Bi,{T)) and neglect the contribution of scattering to the total extinction in the optically 
thin layers [32] • Furthermore, we use a multi-group or opacity-binning approximation \6'i\ I34j to 
account for the dependence of opacity on wavelength: we sort the monochromatic opacities 
into groups (or bins) according to wavelength range and opacity strength, then solve the 
radiative transfer equation for the individual group mean opacities and the integrated group 
source functions [32]. We calibrate the opacity-binning method for each specific choice of stellar 
parameters in order to achieve an as accurate as possible representation of the heating rates. We 
normally do that by finely adjusting the criteria for opacity-bin-membership for each individual 
simulation until the difference between the heating rates computed with opacity-binning and with 
the full monochromatic solution for the average stratification from the 3D model is minimized. 



Table 1. Stellar parameters of some standard stars for which 3D hydrodynamic Stagger-Code 
models have been computed as part of the grid. For all stars, we have assumed a standard solar 
composition with the abundances of metals scaled proportionally to the relevant [Fe/H], with 
a-enhancement of [a/Fe]= -1-0.4 for the metal-poor models. In addition, for IIE1327— 2326 and 
HE0107— 5240, we have also accounted for the peculiar CNO-enhancement of these stars. 



Star 


Tefr/[K] 


log 5/ [cms ^] 


[Fe/H] 


Sun 


5780 


4.44 


+0.0 


Procyon 


6500 


4.00 


+0.0 


HD 140283 


5750 


3.70 


-2.5 


HD84937 


6400 


4.00 


-2.0 


G64-12 


6500 


4.00 


-3.0 


HD122563 


4600 


1.60 


-3.0 


HE1327-2326 


6200 


4.00 


-5.0 


HE0107-5240 


5200 


2.20 


-5.0 



At the present time, a grid of models using a radiative transfer solution with six opacity 
bins is nearing completion. In addition, we have also computed models for some reference 
stars (see Table [TJ using a more refined opacity-binning scheme with twelve opacity bins. We 
plan to eventually extend the twelve-bin opacity binning scheme to the computation of all 
StaggerGrid models. 

To produce the initial 3D snapshot of a simulation for a given set of targeted stellar 
parameters, we take the physical structure from another simulation snapshot previously 
computed for other stellar parameters and scale it appropriately. We do that by looking at 
the ratios of spatial scales and of various important physical quantities such as temperature, 
density, and pressure from ID model envelopes corresponding to the same stellar parameters. 
The scaled model is then allowed to adjust and relax. Spurious p-mode-like oscillations are then 
damped and filtered out to rid the simulation from excess energy caused by the imperfections 
of the scaling procedure and to allow only the natural modes of oscillation to survive. 



Figure 2. Specific entropy per unit mass for a metal-poor red giant simulation: red-orange hues 
indicate high-entropy regions, purple-blue-green hues low-entropy ones. The physical size of the 
box is about 1200x1200 Mm^ horizontally and 450 Mm vertically. The warm inflowing gas from 
the bottom of the simulation domain has constant, high, specific entropy per unit mass. As 
the ascending gas reaches the optically thin layers at the surface, it rapidly cools via radiation 
losses, lowering its entropy. It eventually becomes denser than the surrounding material and 
falls back toward the interior of the star in narrow downdrafts. 

3. Results and applications 

As an illustration of the typical physical structures resulting from 3D simulations, in Fig. [21 we 
show the gas entropy per unit mass in a representative snapshot of a Stagger- Code red-giant 
surface convection simulation. The apparent characteristic surface convective pattern with large 
granules with warm, ascending gas surrounded by an intergranular network of cooler, denser, 
downfl owing material emerges naturally from the simulations without the need for adjustable 
parameters. The first fundamental difference between 3D and ID model stellar atmospheres is 
therefore that the 3D models self-consistently predict the emergence of density and temperature 
inhomogeneities at the stellar surface and their correlation to macroscopic velocity fields. The 
non-linear dependence of the populations of energy levels of atoms and molecules as well as 
of ionization balance and molecular equilibrium on such inhomogeneities ultimately leads to 
appreciable differences between the profiles and strengths of spectral lines generated using 3D 
and ID models, even in those cases where the ID and the average 3D stratifications are not too 
dissimilar from each other [35 j. 

Another difference, especially important in the context of spectral line formation calculations, 
is that ID model atmospheres of late-type stars generally predict a steeper temperature 
stratification as a function of optical depth than 3D simulations at the optical surface and 
in the layers immediately below it. The effects of this are apparent, for instance, in the different 
predictions of centre-to-limb variations (CLV) in the outgoing continuum radiation intensity. In 
the Sun's case, in particular, the predicted CLV can be tested directly against observations, 
showing that the 3D simulations, contrary to ID models, can successfully reproduce such 
variations across the solar spectrum [14j . 

Finally, a third important difference is that 3D stellar surface convection simulations of metal- 
poor late-type stars predict significantly cooler upper-atmospheric temperature stratifications 
than ID models [211 [35] (Fig. [31 left panel). In the 3D hydrodynamic case, the temperature in 
these layers is essentially regulated by two mechanisms: radiative heating due to reabsorption 
of continuum-radiation by spectral lines and adiabatic cooling associated with diverging fiows 
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Figure 3. Left panel, grey shaded area: Temperature distribution as a function of optical depth 
in a representative snapshot from a 3D Stagger-Code simulation of metal-poor red giant; 
blue line: mean temperature stratification; red dashed line: temperature stratification from 
corresponding ID MARCS model. Right: synthetic profiles for the Fe i (neutral iron) spectral 
line at 5110.4 A computed in local thermodynamic equilibrium (LTE) with the 3D (blue line) and 
ID (red dashed line) models of red giant stellar atmospheres shown in the left panel, assuming 
the same Fe abundance. 

above granules. At low metallicities, the contribution of spectral lines to the total opacity in 
the upper atmosphere is relatively weak: the significance of radiative heating in these layers 
is therefore reduced relative to adiabatic cooling and the temperature balance is shifted to 
lower temperatures. Stationary, ID, hydrostatic model atmospheres do not account for such 
adiabatic cooling component associated with gas expansion. The thermal balance in such models 
is regulated purely via heating and cooling by radiation, ultimately resulting in artificially high 
temperatures compared with the 3D case. 

Differences between the temperature stratifications in the upper atmosphere of ID and 3D 
models can amount to ~1000 K in some cases, severely affecting the excitation, ionization, 
and molecular equilibria in those layers. Such temperature differences can have a particularly 
large impact on the strengths of synthetic line profiles from temperature-sensitive species and, 
consequently, on the elemental abundances that can be derived by comparing theoretical and 
observed profiles. As an example. Fig. [3l right panel, shows the predicted profile of a neutral 
Fe line, computed in local thermodynamic equilibrium (LTE) using a 3D model atmosphere 
of a very metal-poor red giant star and its corresponding ID counterpart, assuming the same 
Fe abundance in the two cases. Because of the line's temperature-sensitivity, the resulting 
3D line profile is significantly stronger than the ID one. This also implies that a 3D LTE 
analysis of this line would require in this case a significantly lower Fe abundance than a ID 
LTE analysis to reproduce the strength of a given, observed, line profile. Differences between 
the derived 3D and ID abundances can be of the order of —0.5 dex or even larger in magnitude 
for lines from neutral atoms and molecules [35\ I36| . Based on the experience with previous 
simulations carried out by our group [371 ES] with a predecessor of the Stagger- Code [10] 
and from preliminary analyses with the present StaggerGrid models, we see that the 3D— ID 
LTE abundance differences derived from lines from molecules or neutral atoms typically tend to 
increase in magnitude (i.e., become more negative) when the metallicity decreases, or when the 
effective temperature increases, or, also, when the surface gravity decreases. At the moment, 
however, these results are still based on a limited number of tests restricted to some parts of 
the grid. We are therefore planning to carry out a more systematic comparison of the results 



of abundances determinations with 3D and ID models and study in greater detail the trends of 
the 3D— ID abundance differences with stellar parameters. 

4. Comparison with other models 

The Stagger and CIFIST grids are conceptually similar in terms of basic structure and 
purposes. Both grids also rely on essentially the same opacity sources and use up-to-date, 
realistic equation-of-state packages and input physics. The main differences are in the adopted 
numerical methods (codes), basic resolution, and physical extension of the models (current 
box-in-the-star StaggerGrid models use a higher numerical resolution and generally extend 
down to deeper layers) and in the implementations of radiative transfer and opacity binning. 
The latter, in particular, may be responsible for the apparent differences between the resulting 
temperature stratifications from simulations of metal-poor stars. The StaggerGrid metal- 
poor simulations predict a cooler temperature stratification in the upper-photospheric layers 
compared with analogous CIFIST models computed for the same stellar parameters f38 t [32 l [39] . 
A systematic comparison of the two grids has not been carried out yet, but it is being planned. 
However, it is important to mention that, in spite of the differences between the two grids at low 
metallicity, the StaggerGrid and CIFIST solar surface convection simulations are actually 
in very good agreement with each other, as well as with the current solar simulation by the 
MURaM group [30]. 

5. Summary and outlook 

In light of the results we have presented here, accounting for the differences between 3D and 
ID models is paramount in order to accurately determine elemental abundances and other 
stellar properties from the analysis of stellar spectra. Once the grid will be completed, we 
will therefore compute synthetic spectra for all 3D models and their ID counterparts, covering 
the range from ultraviolet to infrared wavelengths. These spectra will be used for computing 
synthetic colours, deriving stellar parameters and abundances, studying the properties of stellar 
surface convection across the H-R diagram, and for many other applications (see, e.g., Chiavassa 
et al., these proceedings). 

We will also construct and make publicly available average 3D stratifications from the 
StaggerGrid models. The information from the full 3D structures will also be used to provide 
physical constraints to free parameters used in ID analyses such as micro- and macro-turbulence. 
This will facilitate the implementation of the main results from 3D modelling in existing, 
commonly used, ID spectral line formation packages. This would be particularly useful in order 
to systematically and consistently study the combined effects of granulation and departures from 
local thermodynamic equilibrium in late- type stellar atmospheres (see Bergemann et al., these 
proceedings). 

In conclusion, the StaggerGrid will offer a powerful and flexible tool for progressing toward 
precise and accurate analyses of stellar spectra and elemental abundances determinations, and, 
when combined with Gaia-related follow-up surveys, it will provide a significant leap forward in 
our understanding of Galactic chemical evolution. 
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